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1.  Introduction 


This  note  describes  a  method  for  computing  the  answers  to  the  following 
questions  concerning  the  satellites  described  in  the  title: 

(a)  What  is  the  probability  of  having  m  live  satellites  in  orbit 
immediately  after  the  nth  launch? 

(b)  What  is  the  probability  of  having  m  live  satellites  in  orbit  for  the 
first  time  immediately  after  the  nth  launch? 

Question  (b)  arises,  for  example,  when  one  is  planning  the  task  of  putting 
up  an  initial  set  of  m  communication  satellites. 

2.  The  Probability  of  Having  m  Live  Satellites  at  the  nth  Launch* 

Let  the  satellites  be  launched  one  at  a  time,  let  At  be  the  time  between 
consecutive  launches,  let  p  be  the  probability  that  a  launch  will  be  successful, 
and  let  the  lifetimes  of  successfully  launched  satellites  be  distributed 
exponentially  and  have  the  mean  value  t-  The  probability  that  a  successfully 
launched  satellite  will  survive  until  the  scheduled  launch  of  the  next  satellite 
is  then 


c  =  e 


-At/  t 


U) 


Furthermore,  if  there  are  exactly  j  live  satellites  at  one  of  the  launches, 
then  the  probability  that  exactly  i  of  them  will  survive  until  the  next  launch 

is 


(2) 


For  brevity,  we  write  "live  satellites"  for  "live  satellites  in  the  proper 
orbits"  and  "at  the  nth  launch"  for  "immediately  after  the  nth  launch  " 
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J*et  be  the  probability  that  we  will  have  exactly  i  live  satellites 
at  a  launch  if  we  had  exactly  j  live  satellites  at  the  preceding  launch.  When 
we  write 


q  =  1 -p,  d  =  1-c, 


(3) 


the  formula  for 


is 


♦ 


(TT 


i-l^-Ul 


p  • 


<4> 


The  first  term  on  the  right-hand  side  of  (4)  is  the  probability  that  i  of  the 
already  available  j  satellites  survive  and  the  next  launch  fails;  the  second 
term  is  the  probability  that  only  i-1  of  the  already  .available  j  satellites 
survive  but  the  next  launch  succeeds .  Note  the  special  cases 


The  "transition  probabilities"  A 


if  i  *  0. 

(5) 

if  i  *  j+1 , 

(6) 

if  i  > j+1. 

(7) 

form  the  infinite  matrix 
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More  explicitly, 
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*  “  *  *  •  ‘  *  *  *  *  *  Lj 

Let  be  the  probability  of  having  exactly  m  live  satellites  at  the 

kth  launch.  The  probability  distribution  of  live  satellites  at  this  launch  can 
then  be  described  by  the  vector 


(10) 


The  top  component  of  u  is  the  probability  of  having  no  live  satellites  at 
the  kth  launch,  its  second  component  is  the  probability  of  having  exactly 
one  live  satellite  at  that  launch,  and  so  on.  Note  that 


=  0  if  m>k. 
m 


(ID 


The  vector  describing  the  probability  distribution  of  live  satellites  at  the 
(k+l)th  launch  is 


1 

J 

\ 
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(12) 


Our  immediate  problem  ie  to  show  how  the  vector  (12)  can  be  computed  from 
the  vector  (10). 


Consider  the  probability,  ,  of  having  exactly  i  live  satellites  at 

the  (k+l)th  launch.  If  we  knew  that  we  certainly  had  exactly  j  live  satellites 

at  the  kth  launch,  the  probability  of  having  exactly  i  live  satellites  at  the 

(k+l)th  launch  would  be  A...  However,  all  we  know  is  that  the  probability 

of  having  had  exactly  j  live  satellites  at  the  kth  launch  is  u^.  Consequently 

the  possibility  that  we  had  exactly  j  live  satellites  at  the  kth  launch  contri- 
k+1  k 

butes  to  u.  the  term  A..u,  ;  similarly,  the  possibility  that  we  had  j' 

1  J  It 

live  satellites  at  the  kth  launch  contributes  to  u  the  term  »«  an<* 

so  on. 

Accordingly, 


kfl 

u.  = 


k 

c 

j=0 


a  k 

V> 


(13) 


In  view  of  (11),  Equation  (13)  implies  that  the  vector  u^+*  is  the  product  of 
the  matrix  A  and  the  vector  u^;  that  is. 


k+1 

u 


(M) 


although  only  a  finite  portion  of  A  is  in  fact  needed  on  the  right-hand  side 
of  (14),  It  follows  that  u^+^  =  A^u^,  u^+^  =  A^u^,  and  so  on.  As  we  shall 
now  illustrate,  Equation  (14)  provides  a  method  for  answering  the  question 
(a)  stated  in  the  Introduction. 
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An  Example  .  Let  us  use  (14)  to  compute  the  probabilities  pertaining  to 

the  situation  at  the  second  launch.  For  this  purpose  it  is  convenient  to  speak 

of  a  zeroth  launch,  which  is  actually  no  launch  at  all,  so  that  the  probability 

0 

of  having  no  live  satellites  just  after  it  is  unity  and  the  vector  u  has  the 

form  shown  in  (15).  Multiplying  this  u^  by  the  matrix  A,  given  by  (9),  we 

1  12 
get  the  vector  u  in  (15);  and  multiplying  this  u  by  A,  we  get  the  vector  u  . 


It  follows  that  the  probability  of  having  no  live  satellites  at  the  second  launch 
2 

is  q  +dqp,  the  probability  of  having  one  live  satellite  at  that  time  is 
2 

pq+cqp+dp  ,  and  the  probability  of  having  two  live  satellites  at  that  time 
2 

is  cp  .  Note  that  to  obtain  these  results  we  needed  only  a  finite  portion  of  A. 
In  problems  of  this  kind,  the  formal  equation  that  gives  the  probabilities 
pertaining  to  the  situation  immediately  after  the  nth  launch  is 


n  .  n  0 
u  =  A  u 


(16) 


although  actually  only  a  finite  portion  of  A  is  involved. 

2 

Second  Example.  In  the  preceding  example  the  vector  u  was  computed 

on  the  assumption  that  the  outcome  of  the  first  launch  was  unknown.  To 

illustrate  a  different  situation,  let  us  compute  u^  on  the  assumption  that 

the  first  launch  was  successful,  so  that  the  probability  cf  having  one  live 

satellite  at  the  first  launch  was  unity.  The  vector  u*  now  has  the  form 

2 

shown  in  (17),  and  multiplying  it  by  A  we  get  a  new  vector  u  .  Comparing 
this  vector  with  the  u  in  (15),  we  find  that  the  probability  of  having  no 
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live  satellites  at  the  second  launch  is  now  smaller,  and  the  probability  of 


having  two  live  satellites  at  this  time  is  larger. 

3.  The  "First  Time"1  Problem 

The  planning  of  launches  for  a  communication  system  that  uses  m 
satellites  includes  the  consideration  of  two  distinct  tasks: 

(1)  to  put  up  an  initial  set  of  m  live  satellites,  and 

(2)  to  maintain  the  system  by  replacing  the  satellites 
that  fail  after  the  initial  set  has  been  put  up. 

We  shall  deal  here  only  with  the  first  task  and  with  the  probability  (denoted 
below  by  D^)  that  this  task  will  be  completed  in  exactly  n  launches.  In 
other  words,  we  shall  deal  with  the  probability  that  exactly  m  live  satel¬ 
lites  will  become  available  for  the  first  time  at  exactly  the  nth  launch. 

To  illustrate  the  distinction  between  the  probabilities  u°  and  D°  , 

m  m 

we  assume  that  three  live  satellites  are  to  be  put  up  (m  =  3)  and  consider 
the  case  of  six  launches  (n  =  6).  One  possible  sequence  of  the  numbers  of 
live  satellites  available  just  after  the  first  launch,  just  after  the  second 
launch,  and  so  on,  runs  as  follows: 

1,  2,  0,  1,  2,  3  satellites  .  (18) 

Another  possible  sequence  is 


1,  2,  2,  3,  2,  3  satellites  . 


(19) 
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The  probability  takes  into  account  both  of  these  possibilities.  But  in 

■>  ^ 

the  computation  of  the  sequence  (19)  must  be  excluded,  because  it 
implies  that  the  prescribed  task  would  be  completed  in  four  (not  six)  Launches. 


In  the  study  of  the  "first  time"  problem  we  shall  denote  the  matrix 
of  the  transition  probabilities  by  B  rather  than  A,  and  the  probability 
distribution  vectors  by  v's  rather  than  u's. 


The  necessary  and  sufficient  conditions  for  obtaining  m  live  satellites 
for  the  first  time  at  the  nth  launch  are: 


(a)  at  none  of  the  first  n-1  launches  must  the  number  of 
live  satellites  exceed  m- 1 ; 

fb)  at  the  (n-l)th  launch,  m-1  live  satellites  must  be  available; 
and 

(c)  at  the  nth  launch  the  number  of  live  satellites  must  change 
from  m- 1  to  m. 


The  First  n-1  Launches.  We  shall  now  consider  the  first  n-1  launches  and 
shall  take  up  the  nth  launch  later  as  a  separate  problem. 

The  condition  (a)  implies  that  the  v's  whose  superscripts  are  smaller 
than  n  can  have  at  most  m  nonvanishing  components.  Therefore  these 
v's  (unlike  the  u's)  can  be  treated  as  vectors  having  a  finite  number  of 
components.  For  example,  if  n>3,  the  explicit  form  of  v^  is 


3 

v 


(20) 
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The  infinite  matrix  A,  given  by  (9),  includes  transitions  that  violate 
the  condition  (a).  To  eliminate  them,  we  replace  A  by  the  finite  matrix 


where 


■/ 


B  =  /  B 
B 


00 

B01 

•  •  • 

Ddm-1  \ 

(21) 

l 

10 

B11 

•  •  • 

1 

®1,  m- 1 

V 


B 


m 


1,  m- iy 


B„  =  A^,  i»  j  =  0,  1,  .  .  . ,  m- 1  . 
Equation  (14)  must  now  be  replaced  by 


(22) 


k+1  _  k  _ 

v  =  Bv  ,  k^n-2. 


(23) 


In  particular,  if  all  we  know  is  the  outcome  of  the  zeroth  launch,  we  have 

n- 1  _n- 1  0 

v  =  B  v  ,  (24) 

where  v^  is  the  vector  whose  first  component  is  unity  and  whose  other 
components  are  zeros.  But  if  one  knows  how  many  live  satellites  are 
available  immediately  after  some  actual  launch,  he  would  use  the  procedure 
illustrated  in  our  second  example  concerning  the  u’s. 

Note  that  (23)  takes  care  of  all  the  transitions  which  are  consistent 
with  the  condition  (a)  and  ignores  all  the  transitions  which  ?  »*e  not. 

The  nth  Launch.  We  now  turn  to  the  conditions  (b)  and  (c)  on  page  7. 

The  probability  that  (b)  will  be  satisfied  is  the  element  vn  !  of  the  vector 

1  ^  ^ 
v  ,  computed  from  (23)  or  (24).  In  view  of  (6),  the  probability  that  (c) 

m-  i 

will  be  satisfied  is  c  p.  Therefore  the  formula  for  the  probability 
(which  we  have  denoted  by  )  that  exactly  m  live  satellites  will  become 
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available  for  the  first  time  at  exactly  the  nth  launch  ia 

_n  n- 1  m- 1 

Dm  =vm-Jc  P  ,25> 

We  are  indebted  to  D.  H.  Allen,  Pat  Maxwell,  Jr-  ,  and  C.  A.  Moreno 
for  useful  discussions. 


